***************************************************
****----			Table A12	  		   ----****
***************************************************

clear all
cap log close 
set more off

****---- Change directory

global data "/Users/Replication material/Data"

****---- Analysis

use "$data/dataset_replication.dta", clear

global controls Capoluogodiprovincia  population_density_1911 expenditures_1912_PC surplus_tot_1912_PC exp_police_hygiene_1912_PC exp_just_secur_1912_PC  exp_edu_1912_PC exp_public_work_1912_PC share_literate share_literate_male share_literate_female share_male


qui reghdfe gini_log aug18_jan19_sd  $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg11
qui reghdfe gini_log aug18_mar19_sd  $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg12
qui reghdfe gini_log jan18_jun19_all_sd  $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg13
qui reghdfe gini_log jan18_jun19_1_sd  $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg14
qui reghdfe gini_log jan18_jun19_2_sd  $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg15
qui reghdfe gini_log jan18_jun19_3_sd  $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg16
qui reghdfe gini_log jan18_jun19_4_sd  $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg17
qui reghdfe gini_log jan18_jun19_5_sd  $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg18
qui reghdfe gini_log jan18_jun19_6_sd  $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg19

esttab reg11 reg12 reg13 reg14 reg15 reg16 reg17 reg18 reg19, keep(aug18_* jan18_jun19_*) se r2 staraux collabels(, none) ml(,none) cells(b(star fmt (%9.3f)) se(par)) star(* 0.10 ** 0.05 *** 0.01) replace


qui reghdfe gini_log aug18_jan19_d_sd  $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg21
qui reghdfe gini_log aug18_mar19_d_sd  $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg22
qui reghdfe gini_log jan18_jun19_all_d_sd  $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg23
qui reghdfe gini_log  jan18_jun19_1_d_sd $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg24
qui reghdfe gini_log jan18_jun19_2_d_sd  $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg25
qui reghdfe gini_log jan18_jun19_3_d_sd  $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg26
qui reghdfe gini_log jan18_jun19_4_d_sd  $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg27
qui reghdfe gini_log jan18_jun19_5_d_sd  $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg28
qui reghdfe gini_log  jan18_jun19_6_d_sd $controls, cluster(distretto_FE) absorb(i.Zonaaltimetrica i.N_contribuenti_quartile##i.distretto_FE)
est store reg29

esttab reg21 reg22 reg23 reg24 reg25 reg26 reg27 reg28 reg29, keep(aug18_* jan18_jun19_*) se r2 staraux collabels(, none) ml(,none) cells(b(star fmt (%9.3f)) se(par)) star(* 0.10 ** 0.05 *** 0.01) replace
